NAKED SINGULARITIES IN DUST COLLAPSE AS AN EXISTENCE 
PROBLEM FOR O.D.E. AT A SINGULAR POINT 
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Abstract. The final state of the gravitational collapse of a marginally bound dust cloud 
is formulated in terms of an existence problem for the non-linear differential equation gov- 
erning radial null geodesies near the singular point. Rigorous results are proved, covering 
the complete spectrum of the possible initial data. 



1. Introduction 

It is well known that stable, non singular states of superdense matter can exist only if the 
mass of the final object is less than a physical limit, namely the Chandrasekar limit (about 
IAMq) in the case of white dwarfs or the neutron star limit (of the order of 3Mq) in the 
case of neutron stars. For collapsing objects which are unable to radiate away a sufficient 
amount of mass to fall below such limits, no final stable state is available and therefore 
singularities are formed. 

A famous conjecture, first formulated by Roger Penrose [10] and known as the Cosmic 
Censorship conjecture states that a blackhole is always formed in complete gravitational 
collapse of reasonable matter fields. However, if stated without any further mathematical 
assumption, the conjecture is false, since several examples of naked singularities, i.e. solu- 
tions of the Einstein field equations describing singularities not hidden behind an absolute 
event horizon, are known. It is, therefore, of primary importance to understand the math- 
ematical structure of such singularities, with the final aim of reformulating the conjecture 
as a theorem and hopefully prove it. 

Examples of focussing naked singularities in gravitational collapse firstly arose from 
numerical investigations by Eardley [4] and Eardley and Smarr [5], while the first to per- 
form a formal investigation was Christodoulou [1]. In his paper, Christodoulou used a 
fixed point technique to show that the equation of radial null geodesies for a collapsing 
dust ball starting form rest and having a parabolic density profile has a solution meeting 
the singularity in the past, the latter being thus "visible" to nearby observers. Since then, a 
technique has been developed which makes use of L'Hopital theorem to identify existence 
of solutions with finite tangent near the singularity ("root equation" approach, see e.g. [3]). 
In particular, all the possible endstates of the gravitational collapse of spherically symmet- 
ric dust have been obtained in this way [7], as well as the final states of gravitating systems 
of rotating particles known as Einstein clusters [6]. The root equation technique, however, 
proves useful only if the exact explicit solution of the Einstein field equations is known for 
the case at hand. As a consequence, we are still very far from a complete understanding 
of the Censorship problem even in the simple case of spherical symmetry, since very few 
exact solutions are known. In addition, the root equation approach is essentially related to 
the existence of solutions of a specific kind, that is not a-priori guaranteed. 

In this paper, we give a o.d.e. approach to the nature of the singularities in marginally 
bound dust collapse. Using classical techniques we make rigorous, by explicit construc- 
tion, the results obtained previously with the root equation technique. 
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2. Collapsing dust clouds in General Relativity 

A collapsing sphere of dust in General Relativity is described by the famous solution 
which brings the names of Lemaitre, Tolman and Bondi (see e.g. [8]). We concentrate here 
only the case in which the cloud is marginally bound (the velocity is zero at space infinity). 
Using comoving coordinates, the metric is 

ds^ = -dt^ + {R'fdr^ + R^idO^ + sin^ Odtj)'^) 

(we denote by a prime and a dot the partial derivatives with respect to r and t). The function 
R = R{r, t) satisfies the Kepler-Uke equation of motion R? = 2F{r) /R and is therefore 
given by 

(2.1) R{r,t)^r{l-k{r)t)i 

where k{r) = {3/2)^/2F{r)/r^. In the above formulae, F{r) is the initial distribution of 
mass of the cloud (and thus is a positive function). The energy density is given by 

at i = one has e(r, 0) = ^^p- and therefore regularity of the Cauchy data at r = 
implies F « r"^ as r tends to zero. We assume (as usual) the function F{r) to be Taylor- 
expandable near r = (all our results actually hold true also if F is only of class C^). 
Therefore we put 

Fir) = For^ + F^r"+^ + r(r) 

where T{r) is infinitesimal of order greater than or equal to n+4. The physical requirement 
that the density has to be positive and decreasing outwards further imply that Fq is positive 
and Fn is negative. It follows easily that 

(2.2) fc(r) = l-ar"+7(r), 

where 7(r) is infinitesimal of order greater than or equal to n + 1, a is some positive 
constant and, without loss of generality, we have put k{0) = 1. 

The energy density becomes singular whenever R or R' vanish during the evolution. 
Thus, singularities can be of two different kinds: shell crossing, at which i?' vanishes while 
R is non-zero, and shell focusing at which R vanishes. The shell crossing singularities have 
been frequently considered as "weak" although no proof of extensibility is as yet available 
in the hterature. In any case, in most physically interesting situations such singularities do 
not occur, so that we shall concentrate attention here only on the shell focussing case. 

The locus of the zeroes of the function R{r,t) defines the singularity curve ts{r) by 
the relation i?(r, ts(r)) = 0. Due to formula (2.1), wehavets(r) = l/k{r). Physically, 
ts [r) is that comoving time at which the shell of matter labeled by r becomes singular. The 
singularity forming at r = 0, f = fs(0) is called central and, in dust clouds, is the unique 
singularity that can be naked. To see this, we recall that a singularity cannot be naked if 
it occurs after the formation of the apparent horizon. The apparent horizon (th{r), say) is 
the boundary of the region of trapped surfaces and is defined by the equation R{r, th{r) = 
2F{r), that is 

(2.3) t,ir)=U{r)--^kirfr^ 

so that ts{r) > th{r) for any r > 0. 

To analyze the causal structure of the central singularity, observe that, if the singularity 
is visible, at least one outgoing null geodesic must exist, that meets the singularity in the 
past. Such a geodesic will be a solution of 

,2.4) ^ = .(M) 
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where 



(2.5) ip{r,t):=J- 



_ l-k{r)t- ^rk'{r)t 
goo ~ (l-fci)V3 

with initial datum t{0) = ts{0) = 1. For a problem of this kind, in which the initial point 
is singular (the function Lp is not defined at (0, t(0))) no general results of existence/non 
existence are known. As a consequence, in the literature, an approach has been developed 
[3] which makes use of I'Hopital theorem to identify the possible values of the tangent of 
the geodesic curve at the singularity. What turns out is the following: 

• For n = 1 or n = 2 the singularity is naked; 

• For n = 3 the singularity is naked if a > ac where 

_ 2(26+15V3) 

(2.0) ttc = — . 

Therefore, is a "critical parameter": at a = Oc a "phase transition" occurs and 
the endstate of collapse turns from a naked singularity to a blackhole. 

• If n > 3 the singularity is covered. 

This approach, however, strictly depends on the form of the solution of (2.4), that must 
be of the form l + xr" with x constant. Anyway, the root equation can be as well recovered 
following our approach for proving nakedness, where we will look for solutions of the form 
1 + x{r)r", and impose a continuity condition on the unknown function x{r). 



3. Non-existence 

We begin by stating the non-existence result. The argument covers the case n > 4 and 
gives a partial answer in the case n = 3, (the remaining part is given in Section 5). 

3.1. Theorem. If n> 4: the singularity is covered. 

To prove the above statement we need the following: 

3.2. Lemma. There exists r* > such that the apparent horizon th{r) is a subsolution of 
{2 A) for re (0,r*). 

Proof Recall that th{r) = — ^k{r)'^r^, and k{r) = 1 — or", where in last relation 
a > and infinitesimal of order greater than n has been dropped. We must show that 
th{r) < <f{r, th{r)), Vr e (0, r*) with r* sufficiently small. One gets 

(3.1) ^ = -^-1(2^^'^' + 3^'^')' 
and 

(3.2) ^{r, th{r)) = + ^k^r^ + ^kk'r\ 
where the relation 



l-kth= -k'r' 



has been used. It follows 



(3.3) ^ - ^^rMr)) = —kr^ (^^k'r + k 

that is negative for r sufficiently small and positive. □ 



Proof of Theorem 3.1. Let tp{r) the solution of t'{r) = ip{r, t{r)) such that tp{0) = 1. 
By contradiction we suppose the existence of ri > such that tp{ri) < th{ri) and 
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tpi''') ^ th{r), Vr G [0, ri]. We can suppose ri < r*, where r* comes from Lemma 3.2. 
Since tp(0) = t/i(0), one has 

(3.4) < thin) - tpin) = (thin) - thio)) - {(tpin) - t^(o))) = 

where ^ e (0, ri). Using Lemma 3.2 it is t^(^) < <^(^, ift(O)' hence 

(3.5) 4(^)_^(^,f^(^))< ^(^,0) 

Combining 3.4 and (3.5) one gets a contradiction if ^(C, < 0. Now: 

dip, , 2 A;(l - fct) + rfc' (1 - 

(3.6) ^(^.*)-- V 3 ; 



a* ' ' ' 3 (1 - fcf)4/3 

-r^y ^ w lOr small values of r. Showing that ^ 



Hence ^(?', 0) < for small values of r. Showing that ^{r, th{r)) < we have that 



^ i^i S) < 0, since we can observe that the numerator in last term of (3.6) is linear in t 
and the denominator has a fixed sign. Using (3.6) we have 
(3.7) 

d^,^.^, 2^fcV + fcV-|fcfcVt _ f2\-'f8. 1 



where infinitesimal of order greater than n has been dropped out in last quantity, so that 
the sign of the right hand side in (3.7) depends on n, and is strictly negative if n > 4. □ 

The above argument provides only a sufficient condition for nakedness. Indeed, it does 
not exhaust all cases for the singularity to be covered (see Section 5). 

4. Existence 

In this section we estabUsh rigorously existence of naked singularities in the cases n = 
1,2. 

4.1. Theorem. Ifn =1,2 there exists a geodesies of the form 

(4.1) tp{r) = 1 + x(r)r«, r e [0,r*], 

where a = 1 + |n and x{r) is a dijferentiable function in [0, r*] such that x{0) > 0. 

Proof. We will show the existence of a function x{r) E H^'P[0, r*] with p > 1 and r* > 
sufficiently small, and of a parameter a > 1 such that a;(0) > and t = 1 + xr" solves 
equation t'{r) = ip{r, t{r)), where (p{r, t) is given by (2.5): 

^ l-k{r)t{r)-%rk'{r)t{r) „ ,^ 

(4.2) t'{r) = (i_fc,(;))i/3 ■ fc(r) = l-ar" + 7(r). 

We recall that 7(r) is infinitesimal of order greater than n as r ^ 0+. Substituting in (4.2) 
the expressions for t{r) and k{r), and using that n > 1 one gets 



(4.3) rx = r^ 



.1 _ „f n+l-a 



a + fan + rb{r) + xc{r)r°' - d(r).xr""" 
[a + axr'" - xr"-" - S{r){l + xr"))^^^ 



ax, 



where 6(r), c(r) are continuous functions differentiable in r = 0, d{r) = 1 + 7(r), S{r) 
is an infinitesimal differentiable function of order greater than or equal to 1 for r — > 0+, 
and a is a positive parameter to be determined below. We search for x such that r x' is 
infinitesimal for r 0*. Then the right hand side of (4.3) must be infinitesimal for r 
0+. Since the quantity in square brackets is bounded for n = 1,2, it must be |n + 1 — a > 
0. But if the strict inequality held, the Umit of the left hand side would be —ax{0) which 
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by hypothesis is strictly negative. So the only possible situation is a = 1 + from which 
one gets 

(4.4) a;(0) = ai 

using the infinitesimal behaviour of the left hand side of (4.3). 

Having chosen the values of a and the initial condition a;(0), one has actually to show 
the existence of a solution of 



(4.5) rx' 



a [l+^n] +rb{r)+xc{r)r^~^^"- 



d{r)xr^ 



1 + -n ) xG{x,r) 



G-'^{x,r), x{0)=ai, 



where 

(4.6) G(x,r)=a^/3 
= «i/3 



ha;r^+3" -\ \- -xr^+a"^ 

a a 



a 



-5 



xr 3 



S 6 
a a 



We observe that last relation has been written using Taylor expansion of the quantity in 
round bracket in the first row of (4.6), and in view of this the continuous function A{r) has 
been introduced. Using (4.6) in (4.5), and collecting terms with the same power of x the 
differential equation becomes 



(4.7) rx' = G{x,r)-^ 



a{l + -n 



(jl/3 [lJ^_n]j._^_ 



+rb{r) + (e(r) x + f{r)x'^) r^-3"l , x{0) = a^^^ 



where e(r), f{r) are continuous functions, differentiable in r = 0. It is a straightforward 
calculation that e(0) = -1, /(O) > 0. 
With the positions 



,2/3 



/? = a — n = l — -^n, 



y = X — a 

one recovers a differential equation of the form 

(4.8) ry' = A{r,y)y + B{r,y)r^ , y(0) = 0, 

where A{r, y) and i?(r, y) are continuous functions such that >1(0, 0) < and i?(0, 0) > 0. 
Let us now define four constants that bound A and B in a small neighborhood U = 

[0,r*] X [-e,e] of (r, y) = (0, 0): 

Ao<A{r,y)<Ai<Q, < < B{r,y) < B^, {r,y) eU. 
Let us also define the two positive functions 



(4.9) 



^o(r-) 







zi{r) 



respectively solutions of the Cauchy problems 



re [0,r*] 



(4.10) 



z' = I Aoz + Borl^-^ 
z{0) = 



z' = I Aiz + Bir^-\ 



It is readily observed that zo{r) < zi{r) Vr G [0, r*]. Hence, Vn G IN, let y„ denote the 
solution of the ODE in (4.8) with the initial condition y{-^) = yon such that 

Zo{-) < yOn < Zl{-). 
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From comparison theorems in ODE one gets 

(4.11) 0<zo(r) <y„(r)<zi(r), re[l/n,r*], 

and then extending yn to [0,r*] setting y„ = yon in [0, ^] we have that |y„| are equi- 
bounded by K with K constant. Moreover, using the ODE in (4.8) |?/^ | are equibounded 
by K r^^^ which is i^, p > 1. So, up to subsequences, y„ converges uniformly to a func- 
tion y in ff^'P, which is easily shown to be a differentiable solution of (4.8) using the ODE 
in (4.8) and Lebesgue theorem. □ 



5. The critical case 

The analysis so far shows existence of naked singularities if n = 1,2, and non-existence 
if n > 3. When n = 3 a partial answer is contained in Theorem 3.1. Indeed, the key point 
in the proof is the study of the sign of (3.7) for small values of r. In the case fc = 1 — ar^, 
omitting infinitesimal of order greater that 3, direct substitution in (3.7) yields 

t<-«")-(5)"'(|- 

Then we must impose the condition ^ — a > in order to recover the same situation as 
in Theorem 3.1. In other words, we have shown the following 

5.1. Proposition. Ifn = '3 and a < ^ the singularity is covered. 

Sufficient conditions to ensure existence of naked singularity can now be given, with a 
repetition of the argument used in Theorem 4. 1 . In this case one can show the existence of 
a solution of the kind t{r) = 1 + x{r)r^ , with x{r) e H^'P[0, r*], p > 1 and x{0) > 0. 
Since a — n = we must be careful in treating the infinitesimal terms in the differential 
equation (4.3), which now takes the form 



(5.2) rx' = 



3a — d{r)x + rb{r) + xc{r)r'^ 



{a — X + axr^ — S{r){l + xr^)) 



1/3 



- 3a;, 



where b{r), c(r) and d{r) have the same meaning as in (4.3). In order to ensure the infini- 
tesimal behaviour of the right hand side of (5.2), we must then require 

3x(0) = 0. 



(a-x(0))i/3 

Since we want a;(0) > 0, this implies that a must be such that the algebraic equation 
(5.3) 27a;^(a-x)-(3a-a;)^ = 0. 

has real positive roots. It is a simple exercise to check that this is true only if a. < qq or 
a > Uc, where qq = (2/27)(26 + 15v^)~^ while Uc is defined in (2.6). The first case 
however must be excluded since the solution would not live below the apparent horizon 
th{r). Indeed, we know from Proposition 5.1 that the singularity is covered if a < We 
must instead accept the second interval, and the same arguments of Theorem 4. 1 can be 
used with some slight modifications here, in order to ensure the following 

5.2. Proposition. Ifn = 3 and a> ac the singularity is naked. 

What remains to be analyzed is whether naked singularities may exist for a e ( ^ , Oc ) . 
Actually, such solutions represent blackholes, since we can show that the sufficient condi- 
tion of Proposition 5.2 is also necessary in this case. 

5.3. Proposition. Ifn = 3 and the singularity is naked, then a > ac. 
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Proof. Let <p(r) be a solution of the differential equation t'p{r) = Lp{r, tp{r)). We can 
write it in the form tp{r) = 1 + x{r)r'^, although in this case we don't know the behaviour 
of x{r) near the origin r = 0. We just know x continuous, x(r)r'' ^ as r ^ 0+ and, 
since the singularity is naked, tp{r) < th{r) = 1 + {a — ■§f )r^ + o{r'^). Last fact implies 

g 

(5.4) x{r) <a- — +r], 

for r > sufficiently small, where 77 <C 1 is a constant. Then x{r) is bounded from above 
in a right neighborhood of r. But it is also bounded from below. Indeed, t'p{r) > since 
(p{r,tp{r)) > Oforr > small, and then a;(r)r^ is increasing. Thusa;(r)r^ mustapproach 
from above as r — > 0+, and then x{r) must be positive, and therefore bounded. 
Now let us write (5.2) as 

.... / 3a - X + /(r) 

(5.5) r x' = \' - 3x 

(o — a; + g{r, x))'''-^ 

where /(r) — > and, since xr^ — > 0, also g{r, y) ^ per r — > 0+. 
Moreover, let us define 

Q{a,x) = ^. - 3x. 

[a — x)'-''^ 

Recall that (5.4) ensures that (a — x)^^^ > for r small. Then (5.5) may be written as 

(5.6) Q{a,x) — rx' — h{r,x), 

with h{r, x) ^ as r ^ 0+. 

At this point we don't know whether lim^^o x{r) exists or not. If it does, then rx'{r) 
tends to 0. Indeed, from (5.6) we get that rx'{r) tends to Q{a, x{0)); if this quantity was 
not null, then x'{r) would behave like ^ in a right neighborhood of 0, and then x{r) would 
not be bounded, behaving like log r. Then (5(a, x(0)) = 0, which means that a is such that 
x{0) is a positive root of the equation (5.3). 

If limr^o x{r) does not exist, since x is bounded there must exists a sequence (r„ , a;„ = 
x{rn)) with r„ and x'{rn) = as n ^ 00. This shows that {.t„} is such that 
Q{a, Xn) 0. Up to subsequences, {xn} converges to a positive root of (5.3). 

□ 

The fact that we were obUged to divide the analysis on the critical cases into two inter- 
vals of values of a is the mathematical reflection of interesting physical phenomenon [2]. 
In fact, we are using formation of the apparent horizon to obtain non existence. Absence of 
apparent horizon is only a necessary condition for nakedness, and in fact there is a interval 
of values of a for which the singularity is not visible, but the slope of the apparent horizon 
does allow for a geodesic to come out. 
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Note added. After completion of this work we became aware of related independent 
paper by Mena and Nolan [9]. 
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